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SUMMARY 


The Kirchhoff integral formulation, which is often applied in aeroacoustic 
problems involving radiation from surfaces, is evaluated for its effectiveness 
in quantitatively predicting the sound radiated from an oscillating airfoil 
whose chord length is comparable with the acoustic wavelength. A rigid airfoil 
section was oscillated at small amplitude in a medium at rest to produce the 
sound field. Simultaneous amplitude and phase measurements were made of surface- 
pressure and surf ace- velocity distributions and the acoustic free field. Mea- 
sured surface pressure and motion are used in applying the theory, and airfoil 
thickness and contour are taken into account. The result of the investigation 
was that theory overpredicted the sound pressure level by 2 to 5 dB, depending 
on direction. Differences were also noted in the sound-field phase behavior. 


INTRODUCTION 

The determination of sound in a bounded region due to fluctuating loading 
acting over a boundary was first given by Kirchhoff (ref. 1 ) in his integral 
formulation of Huygens' principle. Curie (ref. 2) applied Kirchhoff 's results 
to extend Lighthlll's general theory of aerodynamic sound (ref. 3) to incorpo- 
rate the influence of solid boundaries upon the sound field. Ffowcs Williams 
and Hawkings (ref. generalized Curie's formulation to allow for arbitrary 
convective motion of a surface. 

The aeroacoustic formulations have received wide application. However, in 
most practical problems of interest their use invariably involves assumptions in 
order to reduce the equations to tractable forms. A customary assumption for 
problems in which bodies - airfoils, in particular - encounter subsonic flow or 
undergo subsonic motions is that the linear surface terms of velocity and pres- 
sure dominate the production of the sound field, whereas the contributions from 
Lighthlll's quadrupole term and that of surface viscosity are negligible in com- 
parison (e.g., refs. 5 to 14). With this assvunptlon the aeroacoustic equations 
reduce to Kirchhoff integral formulations. 

The purpose of the present study is to examine experimentally the applica- 
bility of Kirchhoff 's formulation in predicting the radiated sound from an air- 
foil undergoing a particular type of oscillatory motion. The study is motivated 
by the lack of definitive experimental data in the literature establishing 
directly that Kirchhoff 's formulation provides the quantitative relationship 
between the linear surface terms and the acoustic field. Previous experimental 
studies of airfoil noise employing the theory (e.g., refs. 11 to 14) intrinsi- 
cally involved factors which complicate the interpretation of the sound field, 
such as the presence of turbulent and sheared flow and complex distribution of 
surface pressure amplitude and phase. Also, the determination of the surface- 
pressure distribution has normally involved the use of airfoil theory rather 
than direct measurement. Such applications of the Kirchhoff formulation have 



not always produced values that are in good quantitative agreement with the 
noise determined experimentally. 

An experimental approach which provides a controlled method of investigat- 
ing the relationship between surface-velocity and surface-pressure distributions 
and the resultant acoustic field is described herein. The noise source is a 
synimetrical airfoil section which is oscillated at small amplitudes about its 
center of gravity by a resonant vibratory system. To allow definitive compari- 
son, tests were conducted in a stationary medium (no airflow), and only acoustic 
pure tones were generated by the body. For the two particular frequencies exam- 
ined, 301 and 477 Hz, the 45.72-cm (1 8.0-in.) chord airfoil was not small com- 
pared with the acoustic wavelength. These frequencies were chosen because in 
typical aeroacoustic problems the frequencies of interest frequently correspond 
to wavelengths on the order of the chord length (see refs. 15 and 16). 

Simultaneous amplitude and phase measurements were made of the velocity 
and pressure distributions on the airfoil surface and of the acoustic free 
field. The measured surface values of velocity and pressure are used in the 
theory, with airfoil thickness and contour taken into account, to predict the 
phase and amplitude of the acoustic free field. Differences are noted between 
the predicted and experimental values of the sound-field amplitude and phase. 


SYMBOLS 

Values are given in both SI Units and U.S. Customary Units. Measurements 
and calculations were made in U.S. Customary Units. Except when noted, surface 
and sound pressures are presented in decibels (dB = 20 log (Pmis/Pref ^ > where 
Pref = 20 yPa (2 x 10“^ dyn/cm^)). 

b airfoil semichord 

c speed of sound in fluid medium at rest 

Fj^ force per unit area exerted on fluid by solid boundary S in 

Xi-direction 

f test frequency, Hz 

G density of volume source function internal to domain V 

g acceleration due to gravity 

H density of volume source function external to domain V 

1 

£ unit vector in x-|- or y-] -direction 

j unit vector in X 2 ~ or y 2 -direction 

k wave number, w/c 

2 





ko,kg,kx see equation (B5) 


unit vector in X 3 - or 73 -direction 
semispan 

component of Mach vector in direction of r, ? • v/rc 

unit normal vector to surface (drawn outward from surface into fluid) 

specified complex surface pressure amplitude at element m 

absolute pressure 

meain value of pressure in medium 

gage pressure, p - p© 

= I* - y| 

vector between observation point and source point, x - y 
surface of body 


SPP,SPV see equation (7) 


Lighthill’s stress tensor 


observer time 


fluid velocity in the y-frame 


volume region 

velocity of any point fixed in the 5 - frame, 3 y(^,T)/ 3 T 
surface normal velocity, v • n 
observer position vector 


source position vector, C + I v( 5 ,t) dx (the y- and ^-frames are 

•'0 

assumed to coincide at source time t = 0) 
angular displacement of airfoil in oscillatory motion 
sound-field directivity angle (see figs. 4 and B1 ) 

Dirac delta function 

Kroneck'er delta ( 6 ij = 1 for i = j; 6 ij = 0 for i j) 
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e 

— ♦ 

C 

n 

0 

K 

p 

Po 

T 

‘Pa 

‘Pp 

‘Psp 

0) 

V2 

C 3 


unit height of airfoil segment (see fig. B1) 

coordinate system affixed to body 

slope of airfoil segment (see fig. B1) 

angle between n and r (see fig. 1) 

unit chord position of airfoil segment (see fig. B1 ) 

instantaneous value of density 

time-averaged value of density in medium 

viscous stress tensor 

emission time of source 

scalar potentials 

phase angle of surface acceleration at leading edge 

phase angle of surface pressure 

phase angle of sound pressure in acoustic field 


angular frequency, rad/sec 

32 . 


Laplacian operator, 


32 


32 


3y-|2 3y22 


differential operator, i + _L_ j + k 


3yi 3y2 


3y3 


value within is taken at a retarded time T = t - - 

c 


Subscripts: 

b bottom of airfoil segment (see fig. B1) 

i,j components of a vector or summation indices 

m particular element of the airfoil equivalent model (see fig. B1) 

ref constant reference value 

rms 
t 


4 


root-mean- square value 

top of airfoil segment (see fig. B1) 



THEORY 


The aerodynamic sound generated from a region which includes rigid surfaces 
in arbitrary convective motion is given by Ffowcs Williams and Hawkings in ref- 
erence 4. The solution relates the air density fluctuation p - Pq = p(x,t) - Pq 
at position x and time t in terms of space and time derivatives of surface 
and volume integrals over source region y. For the coordinate system shown in 
figure 1, the solution is (from ref. 4) 


4ttc2(p - Pq) = 3 f r dV(c) - ^ f [— rl 

3*i Jy L^l 1 - MrlJ Jg [r| 1 - M^lJ 

+ ^ f r_PoZn_1 dSC?) 

3t Jg Lr|1 - M^lJ 

In the volume integral, T^j is Lighthill's fluid stress tensor related by 
Tij = pu^Uj + [(p - Pq) - c2(p - Pq)J - <?ij 


( 1 ) 


( 2 ) 


The Reynolds stress term pu^uj is usually taken as PoU^uj for low Mach num- 
ber flow. The term 5ij[(P - Pq) ” c^(p - represents the effect of heat 

conduction, which is considered very small for low Mach number flow. The term 
Oij is the viscous stress tensor. 


In the surface integrals, the surface force per unit area F^ is given by 


Fi = n^Cp - Pq) - njPij (3) 

where n^ is the ith component of n, which is the unit surface normal (positive 
from surface into the fluid). The term v„ is the normal surface velocity. 

Equation (1) reduces to that obtained from Curie's analysis (ref. 2) when 
surface convection is not a factor. Equation (1) is regarded as a complete and 
exact description of the sound field. 


The application of equation ( 1 ) to determine the sound produced for a 
particular problem invariably involves assumptions, because the required flow 
parameters normally can only be specified to a limited extent. Consider a slen- 
der body, such as an airfoil section, producing sound due to its own motion 
and/or its reaction to an encounter with a complicated subsonic flow pattern. 

An important assumption for many applications is that the sound is dominated 
by that produced by the net local force on the body and the motion of the body. 
The primary argument used to justify this assumption for particular cases is 
that the quadrupole-type noise source mechanism, represented by the volume inte- 
gral of equation (1), is less efficient than the dipole- and monopole-type mech- 
anisms, as represented in the surface integrals. 

If, in addition to ignoring the contribution of the volume integral, one 
also neglects in equation (1) the viscous components of the surface force, then 
the sound field may be described by the following equation: 
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(4) 


4irp’(x,t) = 4ir[p(5c,t) - PqJ 




PoCVn + P’ 
r|l - Mp 


cos 6 


dS(5) + ( 

• 

j p' cos 9 


5 _ 


dS(5) 


as shown by Farassat (ref. 9). In equation (4), p'(5^,t) is the acoustic pres- 
sure at the observer; within the integrals, p* = p(^,t) - Pq is the surface 
”gage” pressure. 


If the surface is stationary or is undergoing solid-body vibratory motion 
of small amplitude, where |vn| « c, equation (4) may be simplified to 

4TTP’(5?,t) = I L r Po°M [pQ 003 e dS(i^avg) + f Q dS(^avg) (5) 

c 3t r Jg p2 

where ?avg> average value of source coordinate y(^,r), coincides with the 
^-coordinate affixed to the body. Equation (5) is Kirchhoff’s integral formula. 
(See appendix A.) 


MEASUREMENTS 

Description 

Apparatus .- A sketch of the experimental airfoil model is shown in figure 2. 
The rigid ssrmmetrical NACA 0012 aluminum airfoil section (thickness-to-chord ratio 
of 0.12) had a chord of 45.72 cm (18.0 in.) and a span of 30.48 cm (12.0 in.). 

The airfoil was part of a resonant torsional vibration system which also 
included torsional bars with inertia disks. The system was designed to produce 
a sinusoidal oscillation of controllable amplitude about the center of gravity 
of the airfoil section, located 19.05 cm (7.50 in.) from the leading edge of the 
airfoil. As illustrated in figure 2, the vibrational mode used was such that the 
torsional motion at the two disks was 180° out of phase with the rigid-airfoil 
rotational motion. This mode was generated and maintained by input supplied at 
one of the inertia disks by an ordinary shaker unit (15.2 cm (6.0 in.) in diame- 
ter and 15.2 cm long) whose amplitude and frequency were controlled by a fre- 
quency generator and a power supply. 

For the tests two sets of 27.94-cm (11.0-in.) long torsional bars with 
inertia disks having diameters of 30.48 cm (12.0 in.) and 19*38 cm (7.63 in.) 
were used, and each produced the described mode at different tuning frequencies, 
which were respectively 301 Hz and 477 Hz. A photograph of the experimental 
model in the 301-Hz configuration is shown in figure 3. Pretest accelerometer 
surveys verified the mode shape of the sinusoidal motion. 

The tests were conducted in the Langley anechoic noise facility. The chamber 
is approximately 7.6 m (25 ft) long, 7.6 m (25 ft) wide, and 7.0 m (23 ft) high. 
Figure 4 shows the airfoil test assembly suspended from the ceiling in the 
chamber. The center of the airfoil span was positioned 1.83 m (6.0 ft) from 
the acoustic wedges of the floor and about 2.1 m (7 ft) from the wedges of the 
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nearest wall. The shaker was connected to the airfoil assembly by a small- 
diameter steel rod secured at the shaker armature and at an off-center position 
on the inertia disk. In addition to transmitting the sinusoidal input (lateral 
rod displacement on the order of typically 10”3 quj) to the resonant vibratory 
system, the placement of the rod prevented assembly misalignment and sway from 
its vertical position. 

Instrumentation .- A 2-g (0.0643-oz) piezoelectric accelerometer was used 
to determine the rigid-airfoil motion. The accelerometer was mounted in an 
instrument module and positioned 4.92 cm (1.94 in.), measured chordwise, from 
the leading edge. All wiring from the module led out of the test assembly 
through the center of a torsional bar in such a way as to minimize extraneous 
signal noise due to vibration. 

To measure the sound field, fourteen 0.635-cm (0.25-in.) diameter and two 
1.27-cm (0.5-in.) diameter, standard condenser free-field microphones were used. 
The 16 microphones were positioned on a 2.44-m (8-ft) radius microphone boom. 

As shown in the plan view of figure 4, the microphones were placed at 12° incre- 
ments to cover a range of directivity angle 8 of 180°. Care was taken to posi 
tion the microphone diaphragm in the plane perpendicular to the airfoil surface 
at the midspan within 1.3 cm (0.5 in.) and at a 2.210-m (7.25-ft) radius from 
the airfoil center of oscillation within 0.3 cm (1/8 in.). 

The airfoil surface-pressure distribution was determined with piezoresis- 
tive silicon pressure transducers (ref. 17). The devices were chosen because 
of their low acceleration sensitivity, small size, and excellent long-term sensi 
tivity and phase stability. The 0.3175-cm (1/8-in.) diameter diaphragm of the 
transducer has its first resonant frequency at 70 kHz, and with a 10-V dc excita 
tion voltage, the transducer sensitivity is -117 dB(ref. 1 V)/Pa. 

Eight fixed-position transducers and one variable-position transducer mea- 
sured the surface pressure. On one side of the symmetrical airfoil, the fixed 
transducers were flush mounted on the instrument module. These transducers were 
located 3.175 cm (1.25 in.) from the midspan and were positioned at 1.4, 9.7, 
21.2, 31.9, 48.6, 61.8, 78.1, and 88.9 percent of the chord, measured from the 
leading edge (choM length: 45.72 cm (18.0 in.)). 

The spanwise surface-pressure distribution was examined by a transducer 
cemented to a thin steel plate which, in turn, was bonded at each chosen posi- 
tion. As a result of this surface-mounting technique, the transducer diaphragm 
was elevated 0.318 cm (1/8 in.) above the airfoil surface. Although this eleva- 
tion precluded direct surface-pressure measurement , the spanwise distribution of 
surface pressure amplitude and phase was accurately determined by referencing 
data obtained to values measured at flush-mounted transducer locations. The 
movable transducer was positioned along the span in line with the 21. 2-, 78.1-, 
and 88.9-percent-chord flush-mounted transducer locations. These positions were 
at 10.4, 25.0, 33.3, and 41.7 percent of the span, measured from the span center 
line. The 10.4-percent-span positions (3.175 cm (1.25 in.) from the midspan) 
were spanwise symmetrical to the flush-mounted transducer locations. An addi- 
tional position employed was along the line of the flush-mounted transducer at 
locations 97.2 percent of the chord, measured from the leading edge. 
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Experimental Procedure and Calibrations 


Instrument calibrations .- For the surface-pressure transducers the output 
voltage signal e is the sum of the combined signals due to pressure p and 
acceleration a at the diaphragm. The voltage output ratio is 

_e 3 Spp + Saa _ p ^ ^aa 

®ref SpPpgf Pref Spp^gf 

where p and a, and thus e, are complex quantities, with the amplitude and 
phase dependent on experimental conditions. The reference voltage e^gf is 
determined by a chosen reference pressure Pref* value of the ratio Sg/Sp 

has been determined by a two-part laboratory calibration on a random sample of 
four pressure transducers. The pressure sensitivity Sp was found by an acous- 
tic calibration, and the acceleration sensitivity perpendicular to the diaphragm 
plane Sg was determined for the transducers in a vacuum (an atmospheric pres- 
sure of 667 Pa (5 mm Hg)). 

A range of conditions was considered up to a frequency of 2 kHz, up to an 
rms acceleration level of 15g, and up to a sound pressure level of 150 dB. The 
average value of Sg/Sp was 0.506 (mV/g)/(mV/Pa) with a standard deviation of 
0.069. The value was independent of frequency and levels of acceleration and 
pressure. 

In addition to the instrumentation- laboratory calibration for the surface- 
pressure transducers, on-site amplitude and phase calibrations were conducted 
for the accelerometer, microphones, and surface-pressure transducers prior to 
and during the data acquisition period. The calibration instruments included 
a 100-Hz, 1gpuj3 vibration calibrator, a 114-dB sound-level calibrator, and a 
microphone-transducer acoustic calibrator whose calibration frequency is deter- 
mined by a frequency generator. The outputs of these calibrators as a function 
of frequency and sound level were compared with the outputs of several calibra- 
tors of different types, and agreement was found. 

The accelerometer signal was determined to be in phase with the accelera- 
tion. Routine amplitude calibrations incorporated a substitute accelerometer 
of the same type. Cross calibrations were made between the accelerometer and 
its substitute at the beginning and end of the test program. 

The signals from the surface-pressure transducers were found to be in 
phase with the pressure and acceleration variations. The transducers were phase 
matched within 1® at both 301 Hz and 477 Hz. Amplitude calibrations were con- 
ducted prior to each series of tests for 301 Hz and 477 Hz at 120 dB. No mea- 
surable change occurred in the pressure-transducer calibration during the experi- 
mental study. 

Because of the condenser-microphone and preamplifier circuitry, the signals 
from the microphones were inverted with respect to the acoustic pressure at the 
diaphragms. The phase calibrations, which incorporated the same acoustic cali- 
brator used for the surface-pressure transducers, showed that the microphones 
were phase matched within a standard deviation of 2.3® for 301 Hz and within a 
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standard deviation of 6.2° for 477 Hz. The routine 1l4-dB amplitude calibra- 
tions revealed less than a 0.25-dB variation in the calibrations for any one 
microphone. 

Anechoic chamber calibration .- The quality of the acoustic field was exam- 
ined for the two pure-tone frequencies of interest with the use of an apparatus 
(ref. 18) which simulates an ideal acoustic point source. The source emission 
of the apparatus is from a small -diameter opening of a necked-down tube. 

For the tests the airfoil assembly and shaker were removed from the chamber 
and the acoustic driver of the calibrator apparatus was secured to the shaker 
stand. The source was positioned in the center of the vacated airfoil assembly 
region, equidistant to the 16 microphones. The results of the test are shown 
in figure 5, where the directivity of the source as found by the microphones is 
plotted in terms of sound pressure level and phase against microphone position 
(the angle 3). The phase was determined by referencing the microphone signals 
to the acoustic driver input signal. The deviation shown from the mean values 
of sound pressure level and phase indicates the quality of the acoustic field. 

The results show anechoic behavior for the test frequencies within ±1 dB 
sound pressure level and ±12° phase (except for results of one microphone). The 
flagged symbols indicate that the phase was unstable for the microphone at 
3 = 132°. The data for the microphone at 3 = 156° are not included in fig- 
ure 5 because the microphone was replaced during the experimental program. 

Data acquisition and reduction .- For the tests the sinusoidal signals from 
all instruments were appropriately conditioned for recording on a l4-track FM 
analog tape recorder which was operated such that it had a flat frequency 
response up to 10 kHz. 

The primary instruments used to reduce the recorded calibration and experi- 
mental data were two graphic level recorders and a phase meter. To enhance sta- 
bility of the phase measurements, it was found convenient to employ two matched 
dynamic tracking filters. The frequency was tuned by the accelerometer signal 
and the frequency bandwidth used was 20 Hz. The accelerometer signal was used 
to determine the motion of the rigid airfoil section and also as a phase refer- 
ence for all pressure transducer and microphone signals. Determination of the 
surface pressure amplitude and phase incorporated the use of equation (6). 


Results 

Surface- pressure and airfoil-motion measurements .- Four cases are consid- 
ered. For a frequency of 301 Hz, the test airfoil-section motion produced rms 
accelerations of 2.99g and 9.44g at the accelerometer location. These accelera- 
tions dorrespond, respectively, to peak angular displacements a of 82.1 yrad 
and 259.3 yrad. For a frequency of 477 Hz, the rms accelerations were 2.82g 
and 8.91g, which correspond to a values of 30.8 yrad and 97.6 yrad, respec- 
tively. For reference it is noted that the maximiom oscillatory displacement 
was located at the trailing edge, where the rms value was 0.049 mm (1.93 x 10"3 
in.) for the 9.44g case. Thus the maximum oscillatory velocity corresponds to 
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a velocity of 0.131 m/s (0.429 ft/sec), which meets stated requirements for the 
validity of equation (5). 

The results for the eight flush-mounted pressure transducers are given in 
table 1 . Presented are the surface pressures and the phases by which the sinu- 
soidal pressures lead the surface acceleration of the leading-edge region, 
denoted by cpp - These results are also presented in figure 6, where the 

levels of the rms surface pressures are given in decibels. Also shown in fig- 
ure 6 are additional measurement values at the movable transducer position near 
the trailing edge at 97.2 percent of the chord. 

It is seen that for both 301 suid 477 Hz, the surface pressure amplitude 
varies in approximate proportion to, and is approximately in phase with, the 
surface acceleration (positive as directed from the surface into the medium). 

The surface-pressure chordwise distributions are smooth and have minimum val- 
ues and maximum phase change near the center of the sinusoidal oscillation (at 
41.7 percent of the chord). 

The results of the spanwise pressure survey for a 477-Hz case are shown 
in figure 7. As previously discussed, the levels from the slightly elevated 
movable transducer are referenced to values measured at flush-mounted transdu- 
cer locations. The pressure phase is invariant over the span, whereas the pres- 
sure level is essentially uniform over most of the span. The surface pressure 
diminishes in level as the edge of the airfoil section is approached. 

Predictions of the acoust i c field and c omparisons with experimental values. - 
The measured values of surface pressure and velocity may be used in the theory, 
as represented by equation (5), to predict the acoustic field for comparison 
with the experimental values. Sufficient information is known to evaluate equa- 
tion (5) accurately, because the pressure amplitudes on the side of the airfoil 
not containing pressure transducers are identical with values obtained from the 
instrumented side, but the phase of the pressure is shifted by 180°. This is 
true because of airfoil geometric symmetry, motion antisymmetry, and the fact 
that all measured quantities (motion, surface pressure, and sound pressure) were 
sinusoidal with no discernible, harmonics. 

The contour of the actual experimental model is approximated by plane seg- 
ments (see figs. 8(a) and (b)). Because of airfoil symmetry, the corresponding 
segments on each side of the airfoil can be considered simultaneously in the 
calculations. The amplitude and phase of the surface pressure are considered 
uniform over the segments, whereas the chordwise dependence for velocity is 
taken into account for each segment. In using equation (5) to develop the 
model, the orientation and position of the surface with regard to an observer 
in the sound field are accurately defined. 

The airfoil equivalent model as shown in figure 8(b) is comprised of eight 
elements that radiate independently, which is consistent with the integral form 
of equation (5). For an observer in the far field of each individual element, 
equation (5) becomes 
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(7) 


8 

p'(x,t) = SPV + SPP (SPV)n, + (SPP)m 

m=1 

where 

(SPV)j, = 1_ r dS (8) 

41. 3t r 

and 

(SPP)m = J_ I [p*] cos 9 dS (9) 

“ 4ttc ^ J„ r 

where the small contribution from the last integral of equation (5) has been 
excluded (see appendix B) . The symbols SPV and SPP are complex quantities 
and are respectively the sound pressure due to the surface-velocity term and 
that due to the surface-pressure term. 

Details of the evaluation of (SPV)in (SPP)m contained in appen- 

dix B. Far-field approximations for each element m are employed; however, 
the summation procedure rendering p*(x,t) represents a numerical integration 
in which the far field of the airfoil itself is not assumed. 

The prediction of the acoustic field for each experimental case is deter- 
mined by using equation (7) with the values of (SPV)jn and (SPP)^ obtained 
from equations (B8) and (B9), respectively. The experimental results contained 
in table 1 and the airfoil section geometric constants and the physical con- 
stants in table B1 are used to evaluate equations (B8) and (B9). 

In figures 9 and 10 are presented the noise data that were measured simul- 
taneously with the surface pressures and accelerations given in figure 6 (also 
table 1). Presented are the sound pressure levels and corresponding phases by 
which the sinusoidal sound pressures lead the .airfoil-surface acceleration at 
the leading-edge region, denoted by 9gp - as a function of the directivity 
angle 3 . All data are for measurement locations in the plane bisecting the 
midspan of the airfoil section at a radius of 2.210 m (7.25 ft) from the center 
of oscillation. 

Also included in figures 9 and 10 are the corresponding theoretical predic- 
tions determined from equation (7) by using the experimental data from table 1. 
All theoretical levels Include a negative 0.7-dB adjustment to account for the 
spanwise-diminished surface pressure levels near the airfoil-section edges (see 
fig. 7). 


DISCUSSION OF RESULTS 

It is seen that for both 301 Hz (fig. 9) and 477 Hz (fig. 10), the theoreti 
cally predicted acoustic levels are higher than those measured by some 2 to 5 dB 
depending on directivity angle 3. The closest agreements appear near 3 = 90°, 
where it is noted the sound pressure would be expected to be determined for the 
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most part by the total fluctuating airfoil lift, since the retarded times for 
all parts of the airfoil are approximately equal for the measurement position. 

Of equal significance to sound pressure level for evaluation of the 
theoretical-experimental comparison are the phase values. For all cases agree- 
ment may be considered acceptable near 3 = 90°. However for 301 Hz at 3 rang- 
ing from 12° to 132°, the theoretical change in phase is about 30° larger than 
the measured change. For 477 Hz at the same 3 range, the theoretical change 
in phase is about 50° lairger than the measured change. 

Note that the results for the microphone at 3 = 156° were not included 
in the acoustic field calibration (see fig. 5). Therefore, no conclusions may 
be drawn on the basis of the data points for this value of 3. 

In view of the disagreement between the theoretical and experimental 
results, an error analysis was performed (see appendix C). The error analysis 
indicated that the experimental noise data of figures 9 and 10 should be accu- 
rate within ±1 dB in sound pressure level and about ±12° in phase for at least 
14 of the 16 microphones. However, because of the number of independent measur- 
ing systems, only a smaller net error could result in the overall directivity 
determination of sound pressure level and phase. 

The accuracy of the theoretically predicted curves of figures 9 and 10 is 
dependent on that of the surface measurements and the theoretical model. On 
the basis of a cumulative component analysis, with allowances for possible exper- 
imental procedural error, an ample error band to assign these curves should be 
±1.5 dB in level and ±15° in phase. 

From the results of the error analysis, it is concluded that for the case 
considered, the Kirchhoff solution did not accurately render the quantitative 
relationship between the surface-velocity and surface-pressure distributions 
and the resultant acoustic field. 

It is of interest to note that this same trend of overpredicting the sound 
has been found in other quantitative studies, which involved an analogous prob- 
lem of noise from stationary airfoils encountering turbulent flow, as shown in 
references 11 and 12. In contrast, good quantitative results were obtained in 
studies of noise of high-tip-speed rotors (refs. 13 and 14), where the solution 
of Kirchhoff 's problem for surfaces in convective motion was employed using nor- 
mal surface-velocity distributions and calculated steady-loading distributions 
(unsteady loading being ignored). 

This experimental study tests suitability of Kirchhoff 's formulation for a 
particular application; therefore the results cannot be generalized for all aero- 
acoustic applications involving surfaces. The tests were conducted for only two 
frequencies, whose acoustic wavelengths were comparable with the chord length, 
and for only one type of oscillatory airfoil motion. Also, as pointed out in 
reference 1 1 , there are significant conceptual differences in the noise produced 
by surface motion and that produced from a stationary surface opposing fluid 
motion, although the same aeroacoustic theory embodies both. The results do, 
however, serve to question the usefulness of normal surface-pressure measure- 
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merits used in Kirchhoff's formulation to predict noise accurately for arbitrary 
situations involving surface- flow interaction. 

A natural inclination may be to consider that the neglect of viscous effects 
in the Kirchhoff formulation should not give rise to as much disagreement as 
shown in the results. However in another respect the surface-pressure data 
(fig. 6) suggest that viscous effects may not be negligible. The pressure load- 
ings appear to approach constant levels at the edges rather than tending to 
zero. This behavior is consistent with that found in some recent studies of 
unsteady aerodynamics (refs. 19 and 20) demonstrating failure of the classical 
Kutta-Joukowski condition at the edges of airfoils operating at high reduced 
frequency. Such discrepancies, as explained in reference 21, may be due to 
linear viscous effects. If it is indeed true that one must consider airfoil 
theory for unsteady, viscous flow to describe the surface-pressure distribution 
properly, then it may be suggested that one should not ignore viscosity in the 
prediction of the sound field. 


CONCLUSIONS 

For a particular rigid body undergoing vibratory motion, the quantitative 
relationship between surface-velocity and surface-pressure distributions and 
the resultant acoustic free field was examined and compared with that predicted 
from the Kirchhoff integral formulation. 

To allow definitive comparison for pressure levels and phase, the acoustic 
data generated were pure tone and the tests were conducted in a stationary 
medium (no airflow). A symmetrical airfoil section was oscillated at small 
amplitudes at two frequencies: 301 and 477 Hz. These frequencies correspond 

to values of the compressible frequency parameter (wave number times airfoil 
semichord) of 1.3 and 2.0. 

The theory was applied by using measured airfoil surface pressure and 
motion to predict the sound field, which was compared with experimental values. 
It was found that for both test frequencies, 

1 . The applied theory overpredicts the measured sound pressure level by 
2 to 5 dB, depending on directivity angle. 

2. The phase variation, over a large range of directivity angle, is over- 
predicted by the theory. 

3. The closest overall agreement for sound pressure level and phase appears 
near the directivity angle which corresponds to a location perpendicular to the 
airfoil chord line. However, even here differences are significant. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 31, 1977 
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APPENDIX A 


THE KIRCHHOFF FORMULATION 


This appendix gives a summary of the derivation of a solution to the inhomo- 
geneous scalar wave equation. This solution contains the surface integral known 
as Kirchhoff's integral formulation of Huygens' principle, which is shown to be 
equivalent to equation (5), the equation applied in this experimental study. 

Consider the sound inside a given domain as being determined by the sound 
sources inside this domain and by sound which enters from the outside. Mathe- 
matically, one relates the sound field to its sources by integrating an inhomo- 

geneous differential equation. Let represent a scalar potential, let 
G(y,t) be the density of the source function internal to a domain V, and let 

H(y,t) be the density of the source function external to domain V. It is 

assumed that throughout the domain V, the medium containing the source G(j?,t) 
is homogeneous and isotropic and contains no surfaces which would reflect or 
scatter sound waves. Within domain V the scalar function T satisfies the 
equation 

V2T _ J_ ^ = -G(y,t) (A1) 

c2 3^2 


where c is the speed of sound. 

Let domain V be closed and bounded by a regular surface S and let $ 
and ^ be any two scalar potentials whose first and second derivatives are con- 
tinuous throughout V and on S. Given this, Green’s integral is valid: 



(W2$ _ ^ 



dS 


(A2) 


where 3/9n denotes diff erentiation along the normal to S (normal positive 
into V). Let 3? be a fixed observation point within V and 

r = y(x^ - y-|)2 + (x 2 - + (X 3 - y 3 )^ (A3) 

be the distance of the fixed point from a variable point y. (See fig. A1.) 

In equation (A2), ^ is identified with the desired solution of the wave 

equation (A1). The function $ is taken as a function relating the result in 
the field at 5? due to an activity at y . 


The Kirchhoff theory of integration is used to obtain ^ (e.g., refs. 1 

and 22 to 25) from the wave equation (A1), the prescribed conditions of conti- 
nuity (eq. (A2)), and a choice of function To allow the desired solution, 

function $ is conveniently taken as the elementary solution of the homogeneous 
equation 
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V2$ _ 3^ = 0 (A4) 

c2 3t2 

which is 

$ = A 6^t - (A5) 


This is the free- field Green's function (see ref. 7), where A is an arbitrary 
constant, t is time, and 6 is the Dirac delta (or Impulse) function. The 
result is the solution to the inhomogeneous scalar wave equation (A1 ) , also 
known as the Kirchhoff formula. 


y(x,t) 



1[g] dV 



(A6) 


where the symbol [G] = denotes a function taken at retarded time 

T = t - \ c/ 

c 

In equation (A6) the volume integral is a particular solution of the inhomo- 
geneous wave equation which physically represents the contribution of all sources 
G contained within V. The surface integral is a general solution of the homo- 
geneous equation (A4), where the integral extends over S and accounts for all 
sources H located outside V* 


The Kirchhoff formula can be applied for an infinite region V contain- 
ing the observer point i and bounded internally by a closed surface S, The 
restriction formerly introduced - that the surface of integration is a closed 
surface containing 5? - is satisfied by assuming that a portion of the surface 

is displaced to infinity, where its contribution to the integral is zero. (See 
fig. A2.) Such a procedure is valid whenever the Sommerfeld condition is satis- 
fied, that is, when the product Tr is less than some constant, which is always 
the case for the acoustic phenomena considered here. 

In reference 2, Curie employed equation (A6) and the Lighthill formulation 
of aerodynamic sound to determine the influence of solid boundaries upon the 
sound field. 

The surface integral of equation (A6) is what may be described as 
Kirchhoff 's integral formulation of Huygens' principle, as seen in refer- 
ence 1. In equation (A6), if one lets T equal the velocity potential, not- 
ing that the acoustic pressure p’(5?,t) = c^(p - Pq) = PQ3T/3t and velocity 
Vji B ♦ n = then one obtains 

dS + j fp *3 QOS 9 (js (A7) 

Js 


4ir p'(5!,t) = 1 L 
c 5t 


I 


Poc[vn] + [P’3 cos 9 
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where in the integrals, p' is the surface gage pressure and 

cos 6 = -n • = -III 

^ 3n 

which is the cosine of the angle between the surface normal and vector 
? = 5? - (see fig. 1). Equation (A7) is identical with equation (5) and, 
except for the omission of convective effects, is that found by Farassat in 
reference 9. 
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APPENDIX B 


DERIVATION OF THEORETICAL MODEL 


The purpose of this appendix is to show the derivation of the velocity 
term (SPV),jj and the pressure term (SPP)jn for the theoretical model used 
for the prediction of the acoustic field in the experimental section. 


Each of the radiating elements illustrated in figure 8(b) is approximated 
by straight plane segments of which the cross section and its coordinate sys- 
tem are shown in figure B1 . In this figure the origin represents the pivotal 
point of the rigid-airfoil vibration, and the observer at position x is at a 
constant radius ro from the origin and is at an angle 3 with respect to the 
X2~axis. The observer remains in the X2-X3 plane. 

In figure B1 the two segments representing the element m are denoted as 
top and bottom. The coordinate yg is given by 

ys = - ^o\\ 

for the top segment and by 

ys = -b[e - n(5 - Co)] (^1 ^ C 1 C2> (B1) 

f or the bottom segment . Here 
b semi chord of the airfoil 

^ unit position on y2-axis, y2/b 

ri slope of segment m 


e unit height, y3/b, at center of segment m 

The normal n^ for the top segment is 


A 

"t 



1 A 

, ^ k 

Vi + 


whereas for the bottom segment one obtains 


A 

nb 


n j - 1 

\/i + \/i + 


The surface velocity Vjj in equation (8) for element m has linear variation 
in the y2-direction. The prescribed sinusoidal velocity is 

vt = iuxxb exp io)t<^[e - ti(5 - Cq)]3 + 
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for the top and 


Vb = itoab exp iu)t^-[e - n(C ~ 5o^3^ 



for the bottom, where a is the angular displacement. The surface pressure p' 
in equation (9) is assumed to be uniform and sinusoidal over each segment. For 
the top segment of element m, P’ = Pm exp i(Dt, and for the bottom segment, 
p» = -P^ exp icot, because of symmetry of the airfoil and antisymmetry of the 
surface motion. The term P^ is a complex quaintity in order to allow phase 
specification between pressure and surface motion. For the top segment the 
normal velocity Vj, is 


^n,t 


V 


= iaxxb^ exp io)t 


1 - n(C - Co)]l 


n/T 


+ Tl" 


and, correspondingly for the bottom segment. 


Vn,b = ^ • nb = -Vn,t 


For the integrals of equatio ns (8) and (9) the diffe rential surface area 
for the element m is dS = dyi0 + dy2 = dyib^l + vr The coordinate 

yi varies from -5, to %, where 1 is the semispan of the airfoil. The fol- 
lowing is obtained for equation (8): 



Equation (9) becomes 


(SPP)rn = 


_ iba)Pm\/ 


I f^2 


^m\/l + f f 

J_o 4 


'cos exp I ia)\^t 


Mt-?)] 


cos 9b exp ^ico(t - c^)_ 




dC dyi 


(B3) 


The radii r^ and rb are respectively the distance from the differential area 
dS of the top and dS of the bottom segment to the observer at x. The angle 
©b is the angle between the normal at dS for the top segment and the radius 
rt* Correspondingly, ©b is the angle between the normal for the bottom seg- 
ment and rb . 

From the coordinate system illustrated in figure B1 , the center of element 
m at C = Co = (Cl C2^/2 is at a radius rn, from x. The value of rm is 
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given by Tj, = \Jpo^ * r’adius rj, makes an angle Bq, 

with respect to the y 2 ~axis which is related by 


003 6» = !;mij!L_(Cob)iji_£c 

"■ ar.^of 


(B4) 


The distance from a point within element m on the y 2 -axis to x is approxi- 
mated closely by r^ + b(C - cos Bjj,. Therefore, the radius rj^ from a 
point on the top segment to x is 

^t = >^m + - ^o'> &m - - 5o^] si" ^m 

For the bottom segment this radius is 

•'b = I’m + - ^o) ®os 3m + b[e - ti(C - 3 

Rearranging the expressions for r^ and r^ gives 


m 


Ob = rm * ^ - |2(kx - kj) » Sck, - kj) 






where 

k 

^0 

k^ 


wave number, o)/c 
= kbe sin 3 ^ 

= kb cos 3jn 
= kbri sin 3ni 

The cosines of angles 6 ^ and 9^ are approximated as follows: 


cos 9^ s: ^ ^ = 




'm 


cos 05 


6 h = "b . ‘ = - 


‘m 



n 

COS 

3m + 

1 


+ r\^ 

\/l + 


n 

cos 

3m - 

1 


+ 

\f^ + ^ 

for 



COS 

9^;, and 


sin 3 


m 


sin 3 


m 


(B5) 


(B 6 ) 


each point on the segments because ^ « rm> where i is the airfoil semispan 
on the y-i-axis. 

For the integrals of equations (B 2 ) and (B3), only first-order approxima- 
tions are required for r^ and r^ in the denominators. These are determined 
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evaluating the previously obtained values of and at the center of 

each segment, which are defined as ^b,o found to be 


= fm - sin 6 m 
rjj^o ~ *^ni sin 6 m 


(B7) 


To evaluate equations (B2) and (B 3 ), it^is noted^that the values obtained 
for the radii r^ and r^ and the angles and 65 are independent of the 

y-| -coordinate. Also, the only dependence of these values on the variable ^ is 
contained in the last terms of the approximations for r^ and rjj in equa- 
tions (B5). Upon integration and rearrangement, equations (B 2 ) and (B3) become, 
respectively. 


(SPV)m 


-PQ0)2ab2jl(l - n2) gxp [i(cot - krm) 

27T“^ 





^s^*^t,o 


and 


^ 1 ( kv + kg ) n ( ^ Q ) 

Bm 

, . i(kx - kg)n(e + 

1 - 

(kx - ^s^^bjO 

L 1 - J 


^ iUm _ ^^m ^ 
*^t , o ’’b , oj 


(B 8 ) 


rQDD^ _ -Pmkli'bN/l + n2 exp [i(wt - krm)] /Am ©t _ 

21 \~^ — ^'b;^ / 


(B9) 


where 

Am = ?xp (iko) f [-i(kx + kgKCa - ^o)] - e*P + kg)(^1 

kx + kg V 

Bm = [exp [-i(kx - kg )(^2 - ^o)] - exp [-i(kx - kg)(Ci 


- 5o)]} 

- 5o)]) 


Cm = exp exp [-i(kx + kg)(C 2 - ^o>] - exp [-i(kx + kg)(C-| - ^o>]) 

kx + kg J 


and 

Dm = exp (-ikQ)r^^ [-i(kx - kg )(?2 ” ^o^] “ exp [-i(kx - kg)(Ci 
kx - ks I 



These equations, when applied to equation (7), render the prediction of the 
sound field as based on the airfoil equivalent model shown in figure 8 (b). Equa- 
tions (b 8 ) and (B9) are calculated by using experimental results of table 1 and 
the coordinates and slopes of the airfoil equivalent model given in table B1 . 

The prediction ignores the near-field term of equation (5). This term is easily 
accounted for in equation (B9) by multiplying Am by the factor 1 + (1/ikrt,o^ 
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and Bjj by the factor 1 + (l/ikr^ q). However, the contribution of these fac- 
tors is small because the observer at x is considered in the far field of the 
element m. 

A more approximate theoretical model is obtained when one assumes that the 
airfoil thickness is zero ( fig . 8(c)), in which case the values^ of e and ti 
of equation (B1) become zero; this requires ko = kg = 0, cos 9^ = -cos 9b» 

I’t o = ^'b o = •^m* *m = ®m» ^m = ® result, the contribution from 

the (SPV^m term in equation (B8) is zero. 

In addition to the preceding approximation in regard to thickness, if the 
value of kjf multiplied by the width of the element small, then 

the terms A^ and B^ each approach the value ~ total force 

on the element m is 2Pn,(2A)(C2 - Ci)b = F^; substituting this value in equa- 
tion (B9) gives 

(SPP). = cos 8t r A _ rj,\] 

4TTrm L ' ° 

Now if the previously mentioned near- field factors are inserted, equation (BIO) 
becomes 

(Spp) , Fm C03 6t exp [i(o(t - ^ ^ 

rm 

which represents the sound radiated by a simple dipole. 
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ERROR ANALYSIS 


In order to 
and experimental 
experimental and 
comparisons. 


evaluate the accuracy of the comparisons between theoretical 
results in figures 9 and 10, one can examine the individual 
theoretical components and their contribution to the 


Experimental Accuracy 

The expected accuracy of the experimental noise data can be determined 
from the results of the acoustic field test in figure 5. The figure shows a 
scatter error of ±1 dB in sound pressure level and ±12° in phase for 14 of the 
16 independent microphone systems. 

The accelerometer measurements are predicted to be accurate within 6 per- 
cent. This would correspond to an error of ±0.5 dB in level for the calculated 
contribution of the surface-velocity term (eq. (8)) to the noise prediction 
equation. As will be shown in this appendix, the contribution of the surface- 
velocity term is very small compared with that of the surface-pressure term 
(eq. (9)). Therefore, any expected accelerometer measurement error contributes 
very little to the experimental -theoretical comparisons. 

The surface-pressure transducers were found to give very repeatable cali- 
bration and pressure measurement results. In determining a reasonable error 
band for the pressure data, one may look to differences found between transdu- 
cers in regard to the sensitivity ratio S^/Sp, defined in equation (6). On 
the basis of calculations made by using the standard deviation of the calibrated 
values of S^/Sp, the expected error is less than ±0.6 dB in surface pressure 
level and less than ±1.0° in phase for any datum point of figure 6. By using 
this error band, an analysis of conceivable error buildup in the prediction of 
the sound field renders an expected error band of ±1.2 dB in sound pressure 
level and ±10° in phase for the theoretical curves of figures 9 and 10. 


Model Accuracy 

The effect of the degree of theoretical-model refinement was examined and 
some results are shown in figure Cl. In this figure various predictions, based 
on model refinements of figure 8, are given of sound pressure level and phase 
as a function of directivity angle g for two of the experimental cases. No 
experimental noise data are shown. 

It is seen that the predictions using the airfoil equivalent model, used 
for the theoretical curves of figures 9 and 10, are closely approximated by 
using a model with assumed zero airfoil thickness in noise calculations. With 
zero thickness, the value of equation (8) reduces to zero and that of equa- 
tion (9) increases slightly; this renders a total prediction reduction of only 
about 0.5 dB in sound pressure level and a change in phase of less than 2.0°. 
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This shows that for the experimental cases considered in this study, airfoil 
contour refinement may be regarded as unnecessary. 

An even more approximate prediction is obtained when the airfoil is 
replaced by two simple dipoles whose strength is determined by summing measure- 
ments over the surface and whose positions are specified near the centroids of 
the in-phase regions (fig. 8(d)). The calculated results for three -values of 
the angle 8 are shown in figure Cl for the two-dipole approximation and are 
seen to be within 1.5 dB of the values for the more exact models. 

The effect on prediction of increasing the number of radiating elements 
beyond eight was examined. It was found that for surface-pressure distribu- 
tions similar to those of figure 6, additional refinement incorporating more 
elements could result in changes of only about 0.1 dB in predicted noise level. 
Therefore, for the cases considered, eight theoretical model elements are suf- 
ficient for accurate prediction. 

Spanwise compactness is assumed in the model and is allowed because of 
the acoustic- field measuring positions, which were in the plane perpendicular 
to the airfoil surface at the midspan. Analysis has determined that the maxi- 
mum expected error due to this assumption would be a negligible difference of 
10~3 dB in sound pressure level and less than 1° in phase. 

Also, as previously mentioned, there was 0.7 dB subtracted from the theo- 
retical values of figures 9 and 10 to account for spanwise diminishing levels 
near the airfoil edges. This should be an accurate adjustment within 0.1 dB, 
assuming the spanwise level continues to drop near the edge in the manner indi- 
cated in figure 7. 


23 



REFERENCES 


1. Baker, Bevan B. ; and Copson, E. T. : The Mathematical Theory of Huygens' 

Principle. Second ed. Oxford Univ. Press, 1969, pp. 1-41*. 

2. Curie, N. : The Influence of Solid Boundaries Upon Aerodynamic Sound. Proc. 

R. Soc.. (London), ser. A, vol. 231, no. 1187, Sept. 20, 1955, pp. 505-51^. 

3. Lighthill, M. J.: On Sound Generated Aerodynamical ly. I. General Theory. 

Proc. R. Soc. (London), vol. 211, no. 1107, Mar. 20, 1952, pp. 564-587. 

4. Ffowcs Williams, J. E.; and Hawkings, D. L.: Sound Generation by Turbulence 

and Surfaces in Arbitrary Motion. Philos. Trans. R. Soc. London, ser. A, 
vol. 264, no. 1151, May 8, 1969, pp. 321-342. 

5. Amiet, R. K. : Acoustic Radiation Prom an Airfoil in a Turbulent Stream. 

J, Sound & Vib., vol. 41, no. 4, Aug. 1975, pp. 407-420. 

6. Sharland, I. J. ; Sources of Noise in Axial Flow Fans. J. Sound & Vib., 

vol. 1, no. 3, July 1964, pp. 302-322. 

7. Morse, Philip M. ; and Ingard, K. Uno: Theoretical Acoustics. McGraw-Hill 

Book Co., Inc., c.1968, p. 321. 

8. Richards, E. J. ; and Mead, D. J. , eds.: Noise and Acoustic Fatigue in Aero- 

nautics. John Wiley & Sons, Ltd., 1968. 

9. Farassat, F.; Theory of Noise Generation From Moving Bodies With an Applica- 

tion to Helicopter Rotors. NASA TR R-451, 1975. 

10. Farassat, F. ; and Brown, T. J. : Development of a Noncompact Source Theory 

With Applications to Helicopter Rotors. AIAA Paper No. 76-563, 1976. 

11. Clark, P. J. F. ; and Ribner, H. S. : Direct Correlation of Fluctuating Lift 

With Radiated Sound for an Airfoil in Turbulent Flow. J. Acoust. Soc. 
America, vol. 46, no. 3, Sept. 1969, PP. 802-805. 

12. Paterson, Robert W. ; and Amiet, Ray K.: Acoustic Radiation and Surface 

Pressure Characteristics of an Airfoil Due to Incident Turbulence. AIAA 
Paper No. 76-571, July 1976. 

13. Farassat, F.; and Brown, T. J.: A New Capability for Predicting Helicopter 

Rotor and Propeller Noise Including the Effect of Forward Motion. NASA 
TM X-74037, 1977. 

14. Hanson, Donald B. : Near Field Noise of High Tip Speed Propellers in Forward 

Flight. AIAA Paper No. 76-565, July 1976. 

15. Hardin, Jay C. : Airframe Self-Noise - Four Years of Research. Aerodynamic 

Noise, AGARD-LS-80, Jan. 1977, Pp. 6-1 - 6-19. 


24 



16. Commerford , G. L.; and Carta, F. 0.: An Exploratory Investigation of the 

Unsteady Aerodynamic Response of a Two-Dimensional Airfoil at High Reduced 
Frequency. AIAA Paper No. 73-309, Mar. 1973. 

17. IS Pressure Transducers. LQ Series - Thin Line Design. Bull. KPS-LQ3 Ser., 

Kulite Semiconductor Products, Inc. 

18. Norum, Thomas D. : Measured and Calculated Transmission Losses of Sound 

Waves Through a Helium Layer. NASA TN D-7230, 1973. 

19. Satyanarayana , B.; and Davis, Sanford: Experimental Studies of Trailing- 

Edge Conditions on an Oscillating Airfoil at Frequency Parameters of Up 
to One. AIAA Paper No. 77-450, Mar. 1977. 

20. Archibald, F. S.: Unsteady Kutta Condition at High Values of the Reduced 

Frequency Parameter. J. Aircr., vol. 12, no. 6, June 1975, pp. 545-500. 

21. Yates, John E.: On the Singular Role of Viscosity in the Theory of Thin 

Airfoils. Rep. 306 (Contract N00014-76-C-0576) , Aeronaut. Res. Assoc. 
Princeton , Inc . , May 1 977 . 

22. Stratton, Julius Adams: Electromagnetic Theory. McGraw-Hill Book Co., Inc., 

1941, pp. 424-428, 460-464. 

23. Lamb, Horace: Hydrodynamics. Sixth ed. Cambridge Univ. Press, 1932, 

pp. 498-502. 

24. Skudrzyk, Eugen: The Foundations of Acoustics. Springer-Verlag, 1971, 

pp. 489-500. 

25. Malecki, I. (Irena Bellert, transl.): Physical Foundations of Technical 

Acoustics. Pergamon Press, Inc., c.1969, pp. 149-162. 


25 



TABLE 1.- TEST CONDITIONS AND SURFACE PRESSURE RESULTS^ 


Transducer 

number 

Transducer 
location, 
percent chord 

f = 301 Hz; 
a = 82 .'1 yrad 

f = 301 Hz; 
ot = 259.3 Urad 

f = 477 Hz; 
a a 30.8 urad 

f = 477 Hz; 
a = 97.6 prad 

Peak pressure, 

1 ^ml > 

Phase, 

fp - 9a. ^ 
deg 

Peak pressure, 

1 Pml > 

Phase, 

9p - 9a. 
deg 

Peak pressure, 
iPml. Pa 

Phase, 

9p - 9a. 
deg 

Peak pressure, 
1 Pml t 

Phase, 

9p - 9a.‘* 
deg 

1 

1.4 

5.77 

356 

24.47 

356 

7.65 

3 

22.73 

4 

2 

9.7 

5.63 

359 

21.36 

1 355 

6.17 

0 

19.53 

0 

3 

21.2 

5.31 

359 

21.35 

355 

6.64 

6 

20.43 

6 

4 

31.9 

1.58 

: 0 

1 

7.42 

343 

3.09 

22 

9.26 

23 

5 

48.6 

2.91 

161 

11.14 

169 

I 

3.01 

139 

i 

9.62 

140 

6 

61.8 

6.44 

164 

26.27 

I 175 

i 7.34 

I 171 

1 23.53 

171 

T 

78.1 

9.49 

i 164 

' 37.93 

175 

: 11.03 

170 

i 

34.88 

170 

8 

88.9 

10.39 

162 

1 

40.26 

176 

11.96 

1 

171 

38.34 

.1 

173 


^Terminology consistent with that of equation (B9). 

phase by which the surface pressure leads the normal surface acceleration near airfoil leading edge. 



TABLE B1.- CONSTANTS FOR EQUATIONS (B8) AND (B9) 


(a) Physical and geometric constants 


Semichord, b, m (ft) 0.2286 (0.750) 

Semispan, 1, m (ft) 0.1524 (0.500) 

Measurement radius, r^, m (ft) 2.21 (7.25) 

Speed of sound, c, m/s (ft/sec) 344 (1130) 

Mean air density, p^, kg/m3 (lbm/ft3) 1.20 (0.075) 

Acceleration due to gravity, g, m/s^ (ft/sec^) 9.8066 (32.174) 


(b) Coordinates and slopes of airfoil equivalent model 


Element 
number , 
m 

Transducer 
location, 
percent chord 

e 

n 

C2 

?1 

1 

1.4 

0.055 

-0.600 

-0.722 

-0.833 

2 

9.7 

. 100 

-.111 

-.522 

-.722 

3 

21 .2 

.117 

-.050 

-.300 

-.522 

4 

31.9 

.116 

.040 

-.022 

-.300 

5 

48.6 

.106 

.038 

.273 

-.023 

6 

61.8 

.085 

.101 

.569 

.273 

7 

78.1 

.058 

.094 

.843 

.569 

8 

88.9 

.022 

.137 

1.166 

on 

00 
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n 



Figure 1.- Coordinate systems and integration region 
for equations ( 1 ) and ( 4 ) . 
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Figure 2.- Experimental model with illustration of torsional 

mode employed in tests. 
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Figure 3.- Experimental model in the 301 -Hz configuration. 



Figure 4.- Plan view of airfoil test assembly in anechoio chamber 

and side view of assembly. 
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Figure 7.- Results of spanwise pressure survey - surface pressure level and 
phase as a function of spanwise location for f = 477 Hz. Airfoil span: 
30.48 cm (12.0 in.). 



(a) Actual airfoil model. 



(b) Airfoil equivalent model (eight elements). 



(c) Zero- thickness model (eight elements). 
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(d) Two-dipole approximation. 


Figure 8.- Illustrations of actual model and theoretical models used 

in the investigation. 
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Figure 9.- Sound pressure level and phase as a function of directivity angle 
B at a measurement radius of 2.210 m (7.25 ft) for airfoil oscillating 
at f = 301 Hz. 
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Figure 10.- Sound pressure level and phase as a function of directivity angle 
3 at a measurement radius of 2.210 m (7.25 ft) for airfoil oscillating 
at f = 477 Hz. 
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SURFACE S 

Figure A1.- Integration region for Kirchhoff formula. 
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Figure A2.- Integration region for Kirchhoff formula when region V 
is bounded internally by surface S. 
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Figure B1.- Coordinate system used in development of sound field 
formulation of an airfoil equivalent element. 




Figure Cl.- Predicted values obtained by 
of the cases considered in the study 


using various acoustic 
(no experimental noise 


models for two 
data shown). 
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